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ABSTRACT 

A  new  method  for  detection  of  the  hydrogen  resonance  in  a  passive  maser,  has  been  tested  and  experi¬ 
mentally  verified.  This  technique,  which  reduces  the  effect  of  cavity  pulling  on  system  performance,  makes 
use  of  several  amplitude  and  phase  measurements  of  the  combined  transfer  function  associated  with  the 
cavity  and  hydrogen  line.  The  atomic  resonant  frequency,  determined  in  this  way,  has  been  shown  to  be 
essentially  free  from  cavity  pulling.  For  example,  we  have  measured  pulling  factors  40  times  lower  than  those 
measured  using  pure  amplitude  or  phase  techniques  for  the  detection  of  the  hydrogen  resonance.  The  smaller 
cavity  pulling  factor  is  important  in  systems  where  cavity  tuning  errors  are  thought  to  yield  an  intrinsic 
limit  on  overall  clock  performance.  The  technique  can  be  useful  as  a  diagnostic  tool,  or  as  the  hydrogen 
resonance  detection  method  in  an  operational  passive  maser. 


INTRODUCTION 

This  paper  will  show  that  sampling  at  discrete  frequencies  both  the  phase  and  amplitude  of  the  forward 
transfer  function  of  a  passive  hydrogen  maser  can  yield  a  discriminator  characteristic  which  is  essentially 
free  from  cavity  pulling.  This  approach  reduces  the  severe  requirements  placed  on  the  cavity  servo  in  a 
passive  system.  The  nominal  pulling  factor,  when  detecting  by  standard  techniques,  is  proportional  to  the 
ratio  of  the  cavity  Q  to  the  hydrogen  line  Q.  In  our  system,  which  is  a  miniature  passive  hydrogen  maser, 
this  pulling  produces  2  x  10“ 14  fractional  frequency  error  per  hertz  of  cavity  tuning  error.  Therefore,  to 
meet  long  term  stability  goals  of  1  x  10“ 15,  we  require  that  the  cavity,  the  nominal  resonant  frequency  of 
which  is  1420.405Af.Hz,  be  controlled  to  «  0.050  Hz.  Control  of  the  cavity  resonator  to  this  degree  requires 
extremely  stable  and  accurate  measurement  of  the  resonance,  along  with  the  ability  to  tune  the  resonance 
in  a  way  which  does  not  degrade  the  measurement  accuracy.  This  can  be  difficult  to  achieve  and  maintain. 
Utilizing  sampling  of  the  phase  and  amplitude  of  the  forward  transfer  function,  we  have  designed  a  system 
to  lock  to  the  hyperfine  frequency  rather  than  some  specific  characteristic  of  the  overall  transfer  function, 
such  as  the  maximum  of  the  amplitude  response.  This  design  has  yielded  a  reduction  of  cavity  pulling  by 
a  factor  of  40  in  our  present  system.  With  further  development,  we  hope  to  be  able  to  reduce  the  pulling 
factor  to  the  extent  that  the  cavity  servo  will  no  longer  be  critical,  thereby  simplifying  the  implementation 
of  a  passive  system. 


THE  MASER  TRANSFER  FUNCTION 

In  the  paper  by  Viennet,  Audoin  and  Desaintfuscien  [l]  and  also  in  the  paper  by  Lesage  and  Audoin 
[2],  a  useful  model  for  the  steady  state  response  of  a  hydrogen  maser  is  given.  The  model  assumes  that  the 
atomic  medium  obeys  the  Bloch  equations.  The  resonant  microwave  cavity  of  the  model  is  driven  by  both 
the  induced  magnetization  of  the  atomic  hydrogen,  and  by  an  interrogating  field. 

The  scattering  parameters  for  the  maser  are  derived  largely  from  the  the  references  [1]  and  [2].  The 
negative  impedance  concepts  are  well  explained  in  [3],  The  normalized  forward  transfer  function  for  a  two 
port  maser,  driven  by  a  matched  generator,  and  driving  a  matched  load  is  derived  in  [1]  and  [2].  We  have 
used  this  transfer  function  to  fix  the  parameters  of  the  equivalent  circuit  shown  in  figure  1.  The  relevant 
parameters  of  the  equivalent  circuit  are  the  cavity  impedance  Zc  and  the  negative  atomic  impedance  Za.  For 
the  angular  frequency  of  the  probing  signal  u,  near  to  both  the  cavity  (wc)  and  hydrogen  (wo)  resonances, 
one  obtains 
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Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Ze(w)  =  Rc(  1  +  jTcu( w  -  we)) 


a) 

(2) 


It  is  convenient  to  normalize  both  Zc  and  Za  to  Re.  We  define  Zc  =  Zc/R,.  and  Za  =  Za/Rc.  The  coupling 
parameters  Pi  and  p2  are, 


P  i 


Ac 


and 


(3) 


where  Zq  is  the  characteristic  impedance  of  the  external  transmission  lines,  and  Rc  is  the  resistive  part  of 
the  cavity  impedance.  The  ideal  transformer  ratios  are  n2  and  n2.  We  define  the  normalized  generator 
impedance  as  Zg  =  Zg/Zo  and  the  normalized  load  impedance  as  Zt  =  Zi/Zo- 

From  this  equivalent  circuit,  the  complete  set  of  scattering  parameters  for  the  two  port  maser  can  be 
calculated.  The  scattering  parameters  are  found  to  be, 


5  —  Ul1  ~tJ|:)  ~  ifs/(3  +  jT2(a>  -ojq)) 

1  +  jTc(w  -  wc)  -  y^s/(l  +  jT2(w  —  w0)) 

~ _ t'Zl 

21 "  i + m*  -  u>c)  -  Tfs/(i+ -  «0)j 

^  ^22  +  jT,;{w  —  —  Y^s /(I  +  jTi{w  ~  Wo)) 

22  ~  1  +  jTr\w  -  uc)  -  ^/(l  +  jT2{w  -  u0)Y 

s  _ ^12 _ 

12  1  +  jTc(w  -wc)  -  +jT2(w  -  wo))' 


where  the  relation  a/(  1  +  5)  =  Ra/  {Rc  (1  +  Pi  +  Pi))  has  been  used.  The  Vij  are  simply  functions  of  the 
coupling  parameters.  They  are, 


i'll  — 

^22  = 

^21  =  ^12  = 

The  two  critical  frequencies  of  the  system  are  the  ( F 
sition  frequency,  and  the  cavity  resonant  frequency, 


1  —  Pi  +  @2 
1  +  Pi  +  Pi 
1  ~  Pi  +  Pi 
1  +  Pi  +  Pi 

iVPiPi 
I  +  Pi  +  Pi 

=  1  ,  trip  =  0 


) 


(8) 

(9) 

(10) 

(F  =  0  ,  mp  =  0  )  hyperfine  tran- 


Wo  =  the  unpulled  angular  transition  frequency  of  the  atomic  hydrogen 

wc  —  the  angular  resonant  frequency  of  the  microwave  cavity. 

The  relevant  system  time  constants  are, 


Ti  =  longitudinal  time  constant  of  the  atomic  hydrogen 
Ti  —  transverse  time  constant  of  the  atomic  hydrogen 
Tc  =  the  time  constant  of  loaded  cavity  resonator  i.e.  Zg ,  Zi  =  1 
Tc„  =  the  time  constant  of  unloaded  cavity  resonator. 
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Other  relevant  system  parameters  are  defined  as  follows: 


a  =  KQcTiT21  oscillation  threshold  parameter 
a  =  1  at  oscillation  threshold,  and  0  <  a  <  1  for  passive  operation 

S  =  - - —  the  saturation  factor 

1  +  Tj{u>  -  u0 )2 

loaded  Q  of  the  resonant  cavity 
Q  of  the  atomic  line 

coupling  parameter  (11) 


Mo  - 

permeability  of  free  space 

Mb  — 

Bohr  magneton 

Ve  = 

the  cavity  volume 

V  = 

the  cavity  filling  factor 

h  = 

Planck’s  constant  over  2t 

i  = 

net  atomic  hydrogen  flux  in  atoms/second 

b  = 

amplitude  of  the  cavity  magnetic  field. 

js  _ 

k~~wT 


Further  details  regarding  these  quantities  are  given  in  [1],  [2],  [4]  and  [5], 

The  scattering  parameters  are  useful  for  precise  calculation  of  the  effects  of  microwave  measurement 
errors  on  system  performance.  The  principal  errors  of  concern  are  those  due  to  source  and  load  mismatch, 
the  effects  of  lossy  transmission  lines,  residual  reflections  in  the  measurement  path,  limited  measurement 
isolation  or  directivity,  the  effects  of  dispersion  in  the  cavity  resonator,  and  the  effects  of  nearby  cavity 
modes.  All  of  these  imperfections  can  be  analyzed  by  properly  modifying  the  two  port  scattering  parameter 
model  of  the  maser  given  in  equations  (4)  through  (7).  Their  effect  on  long  term  stability  of  the  system  can 
then  be  derived  through  sensitivity  analyses. 

In  a  system  in  which  the  generator  and  load  reflection  coefficients  are  Tg  and  Tj,  respectively,  we  find 
that  the  measured  value  of  S2\  is  given  by 


where  we  have, 


$21, m 


_ S 21 _ _ 

i  -  £nrff  -  522rf  +  rar,(sns22  -  s2lSi2)’ 


(12) 


r  —  Zq  1 

1  q  —  * 

Z„  \ 


and 


Or,  in  terms  of  the  normalized  impedances,  we  obtain 


r  ,  = 

Zi  +  1 


(13) 
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(14) 


5  _  VMl  (  (Za  + 1  )(Z,  +  i)  \ 

21,m  2  \zc  +  Z.+frig  +  fr&i)' 

If  Zg  and  Zi  can  be  considered  to  have  constant  magnitude  and  phase  over  the  frequency  range  of  interest, 
then  Sai,m  is  seen  to  be  inversely  proportioned  to  the  total  loop  impedance  seen  from  inside  the  cavity.  Under 
these  conditions,  we  can  make  the  following  approximation: 


C  ■~— 

^21,m  — 


(15) 


where  the  total  loop  impedance,  Zt  is  given  by, 


ZT  =  Ze  +  Za+0lZa+p2Zi,  (16) 

and  k  is  a  complex  constant. 


COMMON  DETECTION  TECHNIQUES 

Two  means  of  detection  of  the  hyperfine  resonance  are  commonly  considered,  namely,  maximum 
amplitude  detection  and  phase  inflection  detection.  The  case  of  maximum  amplitude  detection  is  fairly 
intuitive.  The  discriminator  characteristic  becomes  the  first  derivative  of  the  magnitude  of  S2i,m  and  we 
can  set  our  crystal  servo  to  lock  to  the  point  where  this  derivative  is  zero.  This  is  equivalent  to  locating  the 
frequency  (u>a)  where, 


B\ZT\ 

dw 


=  0. 

Wj» 


For  the  condition  of  negligible  saturation,  this  criterion  produces  a  pulling  factor  of, 


(17) 


where  the  pulling  factor  Pa  is  defined  as, 


PA-  «« 

Qa(2-a)’ 

(18) 

_  W  -  Wo 

A  —  i 

wc  -  u>0 

(19) 

and  u>a  is  the  detected  frequency.  The  subscripts  A  and  p  differentiate  between  the  different  detection 
techniques. 

If  the  detection  criterion  of  phase  inflection  is  used,  we  must  force  the  servo  to  lock  to  the  frequency 
( u)p )  at  which, 


d2y> 

dui2 


=  0, 


where 

zT(w)  =  |zTHp>K 

This  can  be  shown  to  yield  a  pulling  factor  Pp  that  is  equal  to  the  Q  ratio, 


(20) 


(21) 


As  the  oscillation  parameter  a  approaches  unity  in  (18),  Pa  approaches  Pv.  In  contrast,  Pp  is  indepen¬ 
dent  of  a.  Since  the  spin  exchange  contributions  to  the  overall  relaxation  are  small  at  low  hydrogen  flux  we 
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can  assume  that  Qa  remains  constant  as  a  is  varied  by  changing  the  atomic  flux.  Thus,  when  the  cavity  is 
mistuned  in  a  passive  maser,  the  dependence  of  Pa  on  oc  would  yield  in  the  amplitude  detection  technique 
a  sensitivity  to  flux  variations  not  present  with  the  phase  detection  technique. 

Both  the  amplitude  and  the  phase  detection  techniques  can  be  implemented  by  applying  a  sinusoidally 
frequency  modulated  carrier  to  the  hydrogen  resonance,  and  detecting  the  amplitude  modulation  produced 
on  the  output  signal.  Ideally, there  will  be  no  amplitude  modulation  detected  if  the  carrier  frequency  is 
equal  to  ui0.  The  amplitude  detection  technique  is  implemented  by  modulating  at  frequencies  less  than 
the  bandwidth  of  the  hydrogen  response  whereas  the  phase  technique  is  implemented  by  modulating  at 
frequencies  much  greater  than  this  bandwidth  [6]. 

SWITCHED  FREQUENCY  MODULATION  TECHNIQUES 

We  are  presently  using  switched  frequency  modulation  techniques  to  simulate  various  detection  criteria. 
The  experimental  system  is  diagrammed  in  figure  2.  The  network  analyzer  is  referenced  to  the  slaved 
oscillator  and  is  used  to  produce  signals  at  fixed  frequencies  relative  to  the  oscillator.  Switched  frequency 
modulation  refers  to  the  process  of  stepping  a  CW  signal  through  a  fixed  set  of  discrete  frequencies.  The 
signal  remains  at  a  single  frequency  for  a  fixed  interval  of  time  during  which  some  measurement  is  made. 
Usually  the  amplitude  and/or  phase  of  the  steady  state  response  is  measured  and  then  the  next  frequency 
is  selected  and  the  measurement  is  repeated.  Using  this  technique,  we  build  up  a  set  of  measurements  of 
the  steady  state  transfer  function.  These  measurements  contain  information  about  the  absolute  frequency 
of  the  slaved  oscillator  relative  to  the  hydrogen  resonance  line.  By  using  frequencies  well  separated  from  the 
hydrogen  line,  so  that  only  the  cavity  response  is  sensed,  we  obtain  information  about  the  cavity  resonance 
relative  to  the  crystal  oscillator.  All  this  information  is  used  to  control  both  the  crystal  and  the  cavity, 
closing  both  servo  loops. 

Using  switched  frequency  techniques,  we  have  tried  three  different  schemes  for  detection  of  the  hydrogen 
resonance.  They  are  amplitude,  phase  and  impedance  detection.  The  amplitude  and  phase  detectors  are 
discrete  simulations  of  the  two  detection  schemes  described  in  the  previous  section.  They  are  implemented 
by  sampling  the  transfer  function  and  estimating  either  the  first  derivative  of  the  amplitude  or  the  second 
derivative  of  the  phase  of  S2i,m.  In  both  cases,  the  estimate  of  the  derivative  is  used  as  the  error  signal  in 
the  servo.  It  can  be  shown  that  the  pulling  factor  for  these  two  techniques  is  essentially  the  same  as  their 
continuous  modulation  counterparts. 


IMPEDANCE  DETECTION 

The  impedance  detection  technique  is  implemented  using  switched  frequency  techniques.  As  its  name 
implies,  it  is  based  on  measurements  of  the  total  loop  impedance  Zt .  The  loop  impedance  is  defined  by  (15) 
and  (16).  One  can  obtain  this  impedance  by  inverting  the  measurement  of  S2i,m  and  setting  the  constant 
k  =  1.  In  the  present  model,  we  will  now  assume  that  the  external  terminations  are  well  matched.  That  is, 
Zg  and  Zi  are  real,  frequency  independent  and  equal  to  one.  These  effects  can  be  lumped  into  the  resistive 
part  of  the  cavity  impedance  and  the  external  terminations  assumed  to  produce  cavity  loading  but  no  cavity 
tuning.  In  view  of  these  assumptions,  one  obtains 


Zt  —  Zc  +  Za  +  /?i  +  02i 

*,«,  (22) 

Figure  3  shows  a  plot  of  Zt  in  the  impedance  plane  for  S  —  0,  and  for  the  cavity  mistuned  by  approximately 
half  its  bandwidth.  Angular  frequency  is  the  parameter  running  upwards  along  the  curve.  Note  that  the 
phase  and  magnitude  of  the  impedance  at  any  frequency  are  easily  obtained  by  drawing  the  vector  from  the 
origin  to  the  point  on  the  curve  corresponding  to  that  frequency.  The  measured  magnitude  and  phase  will 
be  the  magnitude  and  phase  of  that  vector. 
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The  curve  is  made  up  of  two  simple  curves,  a  straight  line  running  from  X  =  — oo  to  X  =  -foo,  and 
a  circle  tangent  to  that  line  at  w  =  w0.  The  straight  line  is  due  to  the  cavity  impedance,  and  the  circle  is 
due  to  the  atomic  impedance.  The  relative  pulling  factors  of  the  different  servos  are  easily  seen  from  this 
diagram.  The  frequency  at  which  the  amplitude  servo  will  lock  is  that  point  on  the  circle  (w^)  which  makes 
the  magnitude  of  the  measurement  vector  a  minimum.  The  phase  servo  will  lock  at  the  point  on  the  circle 
(u)p)  at  which  the  measured  phase  is  equal  to  the  phase  of  the  point  where  the  hydrogen  circle  and  cavity 
baseline  touch.  The  pulling  factors  of  the  different  techniques  are  gauged  by  noting  the  relative  distance 
along  the  hydrogen  circle  between  the  point  of  lock  (either  u>a  or  u>p  )  and  u>q.  For  the  parameters  represented 
by  this  particular  impedance  plot,  it  is  apparent  that  the  amplitude  servo  pulling  factor  is  roughly  half  that 
of  the  phase  servo. 

This  graphical  interpretation  of  the  pulling  factor  has  led  us  to  look  for  ways  in  which  we  could  force 
the  crystal  servo  to  lock  directly  to  w0.  This  would  result  in  a  pulling  factor  of  zero  simply  because  there 
would  always  be  zero  distance  on  the  circle  between  the  lock  point  and  u>o  regardless  of  the  cavity  mistuning. 
A  similar  technique  for  locking  to  a>o  is  described  in  [7]. 

With  reference  to  the  impedance  plot  of  figure  4,  the  impedance  detection  technique  is  based  upon  the 
observation  that  a  vector  A  tangent  to  the  hydrogen  circle  will  point  parallel  to  the  straight  line  of  the  cavity 
response  only  if  the  point  of  tangency  is  at  w0.  We  formulate  the  servo  condition  then,  such  that  A  must 
be  parallel  to  B  if  =  w0.  B  is  a  vector  which  points  in  the  direction  of  the  straight  line  representing  the 
cavity  impedance.  The  servo  system  must  be  able  to  measure  the  vectors  A  and  B  and  to  adjust  the  crystal 
frequency  u>x  so  that  A  is  held  parallel  to  B. 

To  estimate  the  vectors  A  and  B,  the  impedance  Zj  is  calculated  at  the  following  four  frequencies: 

Wj  =  u>x  —  6 

wj  =  w*  -  e 

UJy  =  U)x  +  C 

u>4  =  u)x  +  6.  (23) 

Here,  w*  is  the  crystal  probe  frequency  and  e,  6  are  modulation  frequencies  producing  test  frequencies  that 
are  on  and  off  the  line,  respectively.  We  now  form  A  and  B  using  Zx(w)  from  (22).  Let, 


A  =  Zx[^  3)  — 

B  =  Zt(w 4)  ~  Zx(<*i  1).  (24) 

The  condition  that  A  be  parallel  to  B  is  the  same  as  Im(AB *)  =  0,  where  A  and  B  are  taken  to  be  the 
complex  form  of  the  vectors.  We  calculate  that, 


lm[AB* )  = 


8  aTcT^Sefux  —  wp) _ 

(1  +  T2  (wx  —  wo  +  e)2)(l  +  Tj  (wx  —  Wo  ~  *)2) 


BaTeTi6e 

(1  +  {T2t)2)2 


(u)x  -  u„), 


(25) 

(26) 


with  the  approximate  form  valid  for  w*  -w0  <  e.  The  expression  Im(AB‘ )  is  zero  for  uix  =  wo,  independent 
of  the  value  of  the  cavity  frequency.  Therefore,  we  see  that  the  impedance  detection  scheme  should  produce 
no  cavity  pulling. 
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EXPERIMENTAL  VERIFICATION 

We  have  empirically  shown  that  there  is  a  large  reduction  in  cavity  pulling  when  using  the  impedance 
detection  scheme  in  the  experimental  setup  shown  in  figure  2.  The  measurement  of  pulling  factor  was  done 
by  varying  the  cavity  frequency  by  ±2 kHz  about  the  cavity  lock  point  and  measuring  the  resulting  change  in 
output  frequency  when  using  a  given  detection  scheme.  In  order  to  change  detection  schemes  we  simply  hit 
a  softkey  on  the  controlling  computer  (no  hardware  attached  to  the  maser  need  be  changed).  This  eliminates 
the  possibility  that  the  cavity  Q  or  tuning  coefficients  might  change  between  measurements.  We  were  able 
to  measure  a  pulling  factor  of  2  x  10“ 14  per  hertz  of  cavity  tuning  using  amplitude  detection  and  5  x  10~16 
per  hertz  of  cavity  tuning  using  impedance  detection.  This  represents  a  reduction  in  pulling  of  a  factor  of 
40. 

The  residual  pulling  that  we  see  now  is  not  completely  explained.  It  could  be  related  to  saturation  of 
the  hydrogen  line,  delay  in  the  measurement  circuit,  or  dispersion  in  the  sapphire  cavity.  Anything  that 
would  make  the  true  transfer  function  differ  from  the  ideal  of  figure  3  could  cause  some  residual  pulling  and 
is  therefore  suspect.  Efforts  are  now  under  way  to  eliminate  the  major  sources  of  the  residual  pulling  and 
thus  reduce  the  cavity  pulling  by  another  factor  of  ten. 
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f igure  3 . 


Graph  of  Zj  in  the  impedance  plane,  showing 
the  locations  of  CJ0,  CJA  and  6Jp  . 


figure  4.  Impedance  plane  diagram  describing 
the  impedance  detection  method. 
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